Astron. Nachr. / AN 999, No. 88, 789 -|792](2006) / DOI please set DOI! 



On the use of rotational splitting asymmetries to probe the internal ro- 
tation profile of stars. Application to (3 Cephei stars 

J.C. Suarez 1 >*, L. Andrade 2 , M.J. Goupil 3 , and E. Janot-Pacheco 2 

1 Instituto de Astrofisica de Andalucfa (CSIC), Rotonda de la Astronomfa S/N, Granada, 3004, Granada, Spain. 

2 Instituo de Astronomia, Geofisica e Ciencias Atmosfericas da Universidade de Sao Paulo (IAG-USP), Rua do Matao, 
1226, Sao Paulo, Brazil 

3 LESIA, Observatoire de Paris-Meudon, UMR8109, Meudon,France. 



Received 30 May 2005, accepted 1 1 Nov 2005 
Published online later 

Key words stars: Cepheids - stars: rotation - stars: oscillations - stars: fundamental parameters- stars: interiors 

Rotationally-split modes can provide valuable information about the internal rotation profile of stars. This has been used 
for years to infer the internal rotation behavior of the Sun. The present work discusses the potential additional information 
that rotationally splitting asymmetries may provide when studying the internal rotation profile of stars. We present here 
some preliminary results of a method, currently under development, which intends: 1) to understand the variation of 
the rotational splitting asymmetries in terms of physical processes acting on the angular momentum distribution in the 
stellar interior, and 2) how this information can be used to better constrain the internal rotation profile of the stars. The 
accomplishment of these two objectives should allow us to better use asteroseismology as a test-bench of the different 
theories describing the angular momentum distribution and evolution in the stellar interiors. 
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1 Introduction 

The study of internal rotation of stars is one of the main is- 
sues in stellar physics. Rotation is present in almost all the 
stars, and it interacts with other physical processes acting 
in the stellar interior. In particular, understand the transport 
of angular momentum in the interior of stars is crucial to 
correctly and precisely describe the evolution. Turbulence, 
meridional circulation, mixing of elements responsible of 
different fi gradients during evolution, dynamo effects due 
to the presence of magnetic fields, etc., are some of the 
physical phenomena and processes affected by rotation (see 
e.g. Goupil et al. 2005, Goupil 2009, and Goupil & Talon 
2009, for a review on this topic). 

Nowadays, it is possible to probe the internal structure 
of stars thanks to asteroseismology. In what regards rota- 
tion, progress has been made, during the last decades, on 
the knowledge about the rotation-pulsation interaction. Up 
to now, this problem has been tackled using the perturbation 
techniques to compute the stellar oscillations. These meth- 
ods (see e.g. Dziembowski & Goode 1992; Soufi et al. 1998; 
Suarez, Goupil & Morel 2006) are only valid for slow-to- 
moderate rotators (see a review on the effects of rotation 
on stellar p modes by Goupil 2009). For faster rotators, 
the distortion of the stellar structure due to the centrifugal 
is too large and invalids the perturbation techniques. Non- 
perturbative approaches must thus be considered (Lignieres 
et al. 2006; Reese et al. 2006). 
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The present work consider slow-to-moderate rotators, 
i.e. those for which the parameters e = Q/(G Mj R 3 ) 1 / 2 
and = Vt/v n .£ are small, i.e. the stellar structure is not sig- 
nificantly deformed by the centrifugal force, and oscillation 
frequencies are much larger than the angular rotation rate, 
respectively. The failure of the perturbative approach comes 
first for high radial-order frequency modes. We restrict thus 
this study to rotating stars showing oscillations in a rela- 
tively low-order frequency domain (low g- and p modes, 
or mixed modes), like some <5 Scuti stars, 7 Doradus stars, 
solar-like stars and some massive stars like f3 Cephei. 

The heterogeneity of internal structures (changing with 
evolution and different from one star to another) and pro- 
cesses therein, are strong arguments to assume non-uniform 
rotation. It is thus necessary to take into account possible 
variations (in the radial or angular directions) of the an- 
gular momentum transport, and thus in the shape of the 
internal rotation profile. To do so, we use the oscillation 
code FILOU (Suarez 2002; Suarez & Goupil 2008) which 
corrects the oscillation frequencies for up to second-order 
effects of rotation (including near degeneracy effects) in 
presence of radial differential rotation. The study of radial 
differential rotation has also been used for the asteroseis- 
mic studies, e.g. Casas (2006, 2009), Fox Machado (2006), 
Bruntt et al. (2007a, 2007b); Suarez et al. (2005a), Suarez et 
al. (2006b); Suarez et al. (2007), for S Scuti stars, Rodriguez 
et al. (2006a, 2006b), Uytterhoeven et al. (2008), Moya et 
al. (2005), Suarez et al. (2005b), for 7 Doradus stars, or 
Suarez et al. (2009) for j3 Cephei stars. 
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This technique is applied in the present work to analyze 
the asymmetries of the mode splittings due to rotation. In 
particular, we are interested in examining the behavior of 
the rotational splitting and their asymmetries for different 
mode types (g and p), in presence of radial differential ro- 
tation, i.e. to understand physically how variations of the 
internal rotation profile affect the splitting asymmetries. 
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2 Rotational splittings 



Rotational splittings are being used so far, using inversion 
techniques, to determine the internal rotation of stars. The 
most succeeded example of this is the Sun, for which, the 
internal rotation profile is known precisely thanks to helio- 
seismic inversions (Thompson et al. 1996). 

In the framework of the perturbation theory, rotational 
splittings correspond to the first-order in SI (rotational fre- 
quency) correction to the oscillation frequency u) n ,e,m- In 
that case, the Coriolis acceleration (2S1 x v) dominates, and 
the frequency correction term can be written as 



K m (r,0)Sl(r, 9)d6dr 



(1) 



where K m (r, 9) is the rotational splitting kernel, and R the 
stellar radius. We use this expression to define the general- 
ized splitting S m as 



Wl,r. 



2m 



(2) 



For the sake of simplicity, we assume here shellular rotation, 
i.e. Sl(r) = Sl s [l + 170 (r)], so that the splitting becomes 
independent of m, and its kernel can be written as 



(3) 



where £ r and ^ correspond to the radial and horizontal 
components of the eigenfunction, respectively. The density 
of the star (unperturbed) is represented by po, the radial dis- 
tance by r, and A = £(£+1). The denominator of this equa- 
tion is generally known as the mode inertia. When rotation 
is assumed uniform, the rotational splitting simplifies to 



S : 



SI / K(r)dr 
'0 



(4) 



For a shellular rotation, the rotational splitting can be writ- 
ten in the form given by Suarez et al. (2006a), as 



S = mSl s (C L -1- Jo) 



(5) 



where Cl is the Ledoux constant, and Jo is a 77 dependent 
integral (which is null for uniform rotation), defined by 

1 f R 

Jo = t / Vo(r)[yoi+^ZQ 1 -2y iz i-ZQ 1 ] p r 4 dr .(6) 

J Jo 

with Iq being the mode inertia previously defined in Eq.[3] 



Fig. 1 Two types of rotation profile: linear and shellu- 
lar (local conservation of angular momentum). The linear 
profiles were constructed keeping the same rotational fre- 
quency at the stellar surface, but varying it at the core, ap- 
proximately at 100%, 75%, 50%, and 25% of the original 
one vq_ = 2.42. These are compared with a shellular type 
rotation profile with the same rotational frequency at the 
surface. 

3 Rotational splitting asymmetries 

It often observed that rotational splitting are not perfectly 
symmetric (with respect to axisymmetric modes). Mathe- 
matically, asymmetries of splitting can be defined, in its 
generalized form, as 



-4, 



: + ui +m — 2wq . 



(7) 



Doing some simple algebra, it can be shown that splitting 
asymmetries are dependent only upon second-order (in ST) 
terms, particularly on 

2 



A m = m 1 X 2 



w 



(8) 



where expressions for X2 (the Saio form) can be found in 
e.g. Suarez et al. 2006a or Goupil (2009) and references 
therein. Then, it thus possible to construct the kernel for the 
splitting asymmetries, K2 such as 

r-T! 



A, 



m I SI (r)K2(r)r podr 



(9) 



which is implicitly dependent on the rotation profile 770. In 
this work we are attempting to examine K2 in detail. To 
do so, it is necessary to investigate the behavior of A m for 
different mode types and its sensitivity to variations of the 
rotation profile in different zones of the stellar interior. 

4 Probing rotation profiles using splitting 
asymmetries 

The study of the splitting asymmetry kernel is not an easy 
task. Structure and pulsation variables, together with their 
derivatives, make the K2 mathematical expression complex, 
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Fig. 2 Absolute differences AK (top panels) between the kernels calculated for linear rotation profiles and the reference 
one (shellular profile), for the modes gl (left), pi (right). The K subindices correspond the percentage of the core rotation 
rate of the reference model (see Fig. [TJ Bottom panels show the theoretical asymmetries found for the different modes 
and rotation profiles. The horizontal line represents the splitting asymmetry of the gl and pi modes (left and right panels, 
respectively) found for the reference model. In order to identify the models we consider the following representation: 
asterisks for linear profiles, crosses for smoothed shellular rotation profiles, and rhombus for shellular rotation profiles. 



and hamper the interpretation of its variations in terms of 
physics. It is thus rather difficult to invert the rotation pro- 
file by using the complete expression of K2 (see appendix 
of Suarez et al. 2006a). On the other hand, one expects that 
the pulsation terms in the asymmetry kernel are different 
for different mode types (pure p and g modes, and mixed 
modes). It is thus worth studying the impact of variations in 
the rotation profile on the predicted asymmetries for the dif- 
ferent mode types. For the sake of simplicity, we will con- 
sider models typically representative of j3 Cephei stars. 

To do so, we constructed 8 M Q equilibrium models with 
the evolutionary code CESAM (Morel 1997) for two dif- 
ferent rotation profiles: 1) shellular rotation (we used this 
model as the reference), obtained assuming local conserva- 
tion of the angular momentum during the stellar evolution, 
and 2) linear profiles, which drastically varies (with respect 
to reference model) the rotational velocity in all the stel- 
lar interior. We forced all the models to have the same ro- 
tational frequency at the surface (see Fig. [TJ- In this latter 
case, we computed models with rotational frequency of the 



core v — O c /f2 s = v rc f, 0.75 v re f, 0.50 v re f, and 0.25 f re f, 
respectively. 

We first analyzed the impact of the variations in the ro- 
tation profile on the splitting. We used the oscillation code 
FILOU (Suarez 2002, Suarez et al. 2008) to compute the os- 
cillations. Figure|2](top panels) shows the kernel differences 
between the linear rotation profiles models with the refer- 
ence one AKjQ — Kj(r) — K le f(r) predicted for gl and 
pi split modes. As expected, for both type of modes, the 
largest differences are located around the main peaks of the 
corresponding kernels (near the stellar core for gl, and to- 
wards the stellar surface forpl). Even if pi is a mixed mode 
(with energy in both the surface and the core), AKjq are 
larger for gl. This would help to disentangle (and thereby 
identify) low frequency g and p modes, which are typically 
observed in j3 Cephei stars. 

We then calculated, for each triplet, the corresponding 
splitting asymmetry according to Eq.Qand compared them 
with the reference one (Fig. [2] bottom panels). For illustra- 
tion, we included in the figure the results for models com- 
puted with modified shellular-like rotation profiles (simi- 
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lar to the reference model but modifying the rotational fre- 
quency of the core), and smoothed shellular-like rotation 
profiles, used to discard numerical problems with deriva- 
tives near the convective core edge. 

A first glance indicates that the largest relative varia- 
tions from the reference asymmetry corresponds to the gl 
triplet (up to 60% ), and remains around (up to 40%) for the 
pi triplet. As expected, the largest departs from the refer- 
ence asymmetry, A n corresponds to the linear rotation pro- 
file models (asterisks) which overestimate it. In the case of 
smoothed shellular-like models, the results are very close to 
the reference model for both modes with some marginal ex- 
ception. Finally, the asymmetry found for the non-smoothed 
profiles (rhombus) varies significantly for gl, but remain 
very close to the reference asymmetry value (pi triplet). 

Some preliminary conclusions can be extracted: rota- 
tional splitting asymmetries are strongly sensitive to chan- 
ges of the rotation profile, and such a sensitivity is, some- 
how, differential, depending on the type of mode (g, p or 
mixed mode). Physically, the asymmetries are dependent of 
second-order terms in the oscillation frequency that mainly 
account for the distortion caused by the centrifugal force. 
However, we still need to understand how the structure and 
pulsation variables are related to the asymmetry and its vari- 
ations that allow us to construct a simplified kernel for the 
splitting asymmetry (work in progress). Moreover, thanks 
to the very precise data from space missions like CoRoT 
(Baglin et al. 2003), Kepler (Gilliland et al. 2010), it is pos- 
sible to measure with unprecedent precision the rotational 
splittings and their asymmetries. This, together with the in- 
crease of the number of detected modes (see e.g. Garcia 
Hernandez et al. 2009, Poretti et al. 2009), will definitely 
help for the present study. 
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